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or  any  other  person  or  corporation,  or  conveying  any  rights  or  permission  to 
manufacture  use,  or  sell  any  patented  invention  that  may  in  any  way  be  related 
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FOREWORD 


The  work  described  in  this  report  was  performed  during  fiscal  year  1 976 . It 
was  funded  jointly  by  the  Air  Force  Rocket  Propulsion  Laboratory  under  Job 
Order  No.  305910MB  and  by  Strategic  Systems  Project  Office  under  Naval 
Weapons  Center  Job  Order  No.  14571110.  This  effort  was  part  of  a task  per- 
formed jointly  -by  AFRPL  and  NWC  to  evaluate  structural  computer  codes  for 
analysis  of  rocket  nozzles. 

The  work  was  documented  in  two  volumes.  Volume  I contains  a final  report  and 
summary  of  computer  codes  evaluated  in  this  study;  this  work  was  performed 
mainly  at  NWC.  Volume  II  contains  results  from  solving  several  sample 
problems  using  selected  computer  codes;  this  work  was  performed  mainly  at 
AFRPL. 

This  report  has  been  reviewed  by  the  Information  Office/DOZ  and  is  releasable 
to  the  National  Technical  Information  Service  (NTIS).  At  NTIS  it  will  be  avail- 
able to  the  general  public,  including  foreign  nations.  This  report  is  unclassified 
and  suitable  for  public  release. 
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SECTION  I 


INTRODUCTION 

The  art  of  performing  nozzle  structural  analysis  for  solid  rockets  has  become 
increasingly  more  difficult  in  the  past  several  years.  The  difficulty  has  occurred 
due  to  the  introduction  of  exotic  carbon  materials  to  withstand  increasing  hostile 
nozzle  environments.  These  newer  materials  have  generated  four  problem 
areas  for  the  nozzle  community;  a lack  of  accurate  material  property  data,  a lack 
of  adequate  failure  criteria,  difficulty  with  consistent  material  reproducibility,  and 
a lack  of  adequate  models  in  current  computer  codes  to  analyze  these  materials. 

Substantial  money  and  effort  has  been  and  continues  to  be  spent  by  the  Naval 
Weapons  Center  (NWC)  and  the  Air  Force  Rocket  Propulsion  Laboratory  (AFRPL) 
addressing  the  problems  of  material  properties,  failure  criteria,  and  repro- 
ducibility. However,  little  formal  effort  has  been  directed  toward  eliminating 
the  problems  with  state-of-the-art  computer  codes  relative  to  the  newer 
materials.  This  program  is  a first  step  toward  the  solution  of  this  problem. 

More  than  anything  else,  this  program  confirmed  two  accepted  facts; 

1.  There  is  not  yet  available  a code  which  meets  all  of  the  needs  of  the  nozzle 
community,  2.  The  present  codes  perform  equally  well  for  linear-elastic 
analysis  which,  unfortunately,  is  the  only  type  of  analysis  this  program  had  time 
to  investigate.  One  significant  item  was  discovered.  The  core  requirements 
and  runtimes  for  a general  3-D  code  running  in  the  axisy mmet ric  mode,  in  this  case 
NEPSAP,  were  not  significantly  greater  than  SAASIII  for  analyzing  2-D  axisymme- 
tric  problems.  This  means  that  the  cost  of  using  a 3-D  code  is  not  prohibitive 
for  solving  geometrically  simpler  problems. 

Volume  I of  this  report  presents  a summary  of  nonlinear  analysis,  a dis- 
cussion of  the  features  now  needed  for  stress  analysis,  and  summaries  of  the 
codes  which  were  surveyed  and  evaluated.  This  volume  will  present  the  instal- 
lation of  the  4 codes  which  were  selected  for  running,  the  analyses  which  were 
performed  with  these  codes,  the  data  and  results  of  the  analyses,  and  conclusions. 
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SECTION  II 


INSTALLATION  AND  MODIFICATIONS 

Four  codes  were  chosen  and  installed  on  the  AFRPL  computer  for  further 
evaluation.  All  four  codes  required  some  modification  for  the  purposes  of  this 
evaluation.  Most  of  the  modifications  were  to  reduce  the  core  requirements  so 
that  they  would  fit  within  the  somewhat  small  allowable  core  size  of  the  AFRPL 
computer  (140, OOOg*  locations  or  approximately  49,000  decimal  locations). 

The  four  codes  installed  for  this  program  were  SAAS  III,  TEXGAP,  NONSAP, 
and  NEPSAP.  Each  of  the  codes  will  be  discussed  separately  below  concerning 
version  of  code  used,  modifications,  installation  problems  and  core  required 
for  execution. 

The  SAAS  III  code  which  was  used  for  this  program  was  already  operational 
at  the  AFRPL.  The  code  was  obtained  from  Aerospace  Corporation  in  April  of 
1973.  Since  that  time  the  code  had  undergone  several  modifications  to  reduce 
the  required  core  for  execution.  Another  modification  changed  the  method  for 
interpolating  the  nodal  temperatures.  This  modification  changed  the  scheme  by 
which  the  input  temperature  points  were  chosen  for  interpolation  so  that  problems 
encountered  due  to  high  thermal  gradients  were  eliminated.  Finally,  specifically 
for  this  project,  a time  log  was  added  to  output  the  execution  times  for  various 
segments  of  the  code.  In  this  configuration  SAAS  III  loads  in  125, OOOg  locations. 

The  TEXGAP  code  was  also  operational  before  the  start  of  the  program. 

This  code  was  developed  by  Dr.  Eric  Becker  and  Dr.  Robert  Dunham  at  the 
University  of  Texas  at  Austin  under  an  AFRPL  contract.  The  version  used  for 
this  project  incorporated  modifications  through  Change  4 as  released  by  Dr.  Becker 
and  Dr.  Dunham.  Only  two  additional  minor  changes  were  made  for  this  program. 
The  first  allowed  a two  card  input  for  orthotropic  material  properties  and  the 
second  was  the  installation  of  a time  log  to  output  execution  times  for  various 
segments  of  the  program.  TEXGAP  loads  in  106,000g  locations. 

NONSAP  was  obtained  from  the  University  of  California  at  Berkley  in  June 
of  1975.  The  only  change  made  in  this  code  was  the  removal  of  the  RETURN 
statements  in  the  secondary  level  overlay  programs.  NONSAP  loads  in  131,000^ 
locations. 

*Octal  numbers  are  designated  by  the  subscripted  8. 
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NEPSAP  was  obtained  from  the  Aerospace  Structures  Information  and 
Analysis  Center  (ASLAC)  at  Wright-Patterson  AFB,  Ohio  through  the  Naval 
Weapons  Center.  This  version  was  set  up  for  the  CDC  system, but  it  also 
contained  all  the  required  statements  for  the  UNIVAC  system.  The  UNIVAC 
statements  had  been  essentially  removed  by  making  them  comment  statements. 

NEPSAP  required  substantial  work  to  make  it  operational.  The  biggest 
problem  was  to  reduce  its  size  so  the  code  would  load  into  the  AFRPL  system. 
To  accomplish  this,  all  of  the  I/O  buffer  lengths  were  reduced  to  512  locations, 
and  the  blank  common  length  was  reduced  from  15000  to  12900.  This  allowed 
the  code  to  load  in  127,000g  locations.  Several  other  modifications  were  also 
made.  All  of  overlay  3 was  removed.  This  overlay  set  up  the  instructions  for 
the  post  processing  plot  program  which  the  AFRPL  does  not  have.  The  time 
log  subroutines  had  to  be  rewritten  to  make  them  work  on  the  AFRPL  CDC  6400 
Computer.  In  addition  several  subroutines  still  contained  active  statements  for 
the  UNIVAC  system  plus  many  other  statements  not  required  by  the  CDC  system 
which  had  to  be  removed. 

During  the  course  of  the  program  it  was  discovered  that  NEPSAP  did  not 
have  an  orthotropic  material  model  for  the  2-D  axisymmetric  element.  At  that 
time  a subroutine  from  SAP  IV  was  added  to  NEPSAP  which  partially  fixed  the 
problem.  The  fix  was  partial  because  it  affected  only  the  displacements  and 
strains.  The  amount  of  work  needed  to  make  the  orthotropic  model  work  for 
the  stresses,  due  to  the  nature  of  the  solution  process,  was  determined  to  be 
outside  the  scope  of  this  project. 
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SECTION  III 


ANALYSES 

PROBLEM  DESCRIPTIONS 

The  problems  chosen  for  analysis  with  each  of  the  codes  installed  were 
presented  in  Volume  I of  this  report.  They  are  repeated  here  for  convenience. 

Three  of  the  four  problems  are  problems  which  were  sent  to  all  solid  pro- 
pellant contractors  by  the  JANNAF  Structures  and  Mechanical  Behavior  Working 
Group  and  the  Operational  Serviceability  Working  Group  to  get  a comparison  of 
the  finite  element  capabilities  of  the  various  contractors.  These  problems  were 
chosen  for  three  reasons:  (1)  they  are  practical  problems,  (2)  the  results 

published  in  the  Thiokol  report*  to  the  JANNAF  committees  provides  an  additional 
data  point  for  comparison,  (3)  the  JANNAF  committees  and  contractors  are 
familiar  with  these  problems  which  will  aid  in  their  assessment  of  the  results 
presented  in  this  report.  Results  from  the  Thiokol  analyses  appear  in  this  report 
in  two  cases  for  problem  number  3. 

For  this  program,  two  of  the  JANNAF  problems  were  run  for  two  values  of 
Poisson's  Ratio,  0.4999  and  0.  30.  This  test  was  done  so  that  the  nonrefo rmulated 
finite  elements  used  in  SAAS  III,  NONSAP,  and  NEPSAP  would  not  be  penalized 
for  using  values  of  Poisson's  Ratios  close  to  0.50.  Also  for  problem  number  2, 
sphere  under  internal  pressure,  the  radius  was  not  allowed  to  go  to  zero,  because 
a zero  radius  causes  severe  problems  for  finite  element  solutions. 

Problem  number  1 is  a plane-stress  cylinder  subjected  to  internal  pressure 
only.  The  grid  and  boundary  conditions  which  were  used  for  all  codes  are  shown 
in  Figure  1 . 

The  other  input  parameters  are: 

Internal  Pressure  = 100  psi 
E = 1000  psi 
V = .4999 

Problem  number  2 is  an  internally  pressurized  sphere.  The  grid  is  basically 
the  same  as  that  used  for  problem  number  1.  The  two  minor  differences  are 

Anderson,  G.  P.  , and  M.  D.  Black,  "Finite  Element  Stress  Analysis  Check 
Problems,"  TWR-7808,  Thiokol  Chemical  Corporation,  Wasatch  Division,  1974. 
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that  the  radius  is  not  allowed  to  go  to  zero,  and  the  radial  displacement  boundary 
condition  on  the  Z axis  was  removed.  This  difference  is  shown  in  Figure  2.  The 
remaining  input  parameters  are: 

Internal  Pressure  = 100  psi 

E = 1 000  psi 

V = .4999  and  . 3 

The  third  problem  is  essentially  a propellant  grain  subjected  to  thermal 
loads  due  to  cooling.  The  grid  and  boundary  conditions  are  shown  in  Figure  3. 
The  other  input  data  are: 

E = 1 000  psi 
adT  = - . 01 

V = .4999  and  . 3 

Problem  number  4 was  taken  from  a TRW  report1'  for  Atlantic  Research 
Corp.  This  problem  is  a cool-down  analysis  of  a SiC/PG  coated  throat  insert 
(no  backup  ring).  This  problem  was  chosen  because  it  is  a true  nozzle  problem 
and  employs  one  of  the  state  of  the  art  materials  alluded  to  previously.  Also 
the  material  properties  of  both  the  ATJ  substrate  and  the  SiC/PG  materials  are 
reasonably  well  defined  to  permit  analysis. 

The  grid  and  boundary  conditions  used  are  shown  in  Figure  4.  The  material 
properties  for  ATJ  and  SiC/PG  are  given  in  Tables  1 and  2.  These  property 
values  were  taken  from  Tables  I and  II  of  the  TRW  report. 

RESULTS 

The  results  of  the  four  analyses  are  presented  in  this  section.  Due  to  the 
large  amount  of  data  that  was  generated,  much  of  it  has  been  omitted  from  this 
report  for  the  sake  of  brevity.  The  data  is  recorded  in  two  forms,  tabular  and 
plots,  and  is  located  in  Appendix  A.  The  plots  provide  a clear  visual  interpreta- 
tion of  the  results  and  show  trends  and  comparitive  accuracies.  The  tabular  data 


’''King,  K.  R.  , "Nonlinear  Stress  Analysis  of  Codeposited  Silicon  Carbide- 
Pyrolytic  Graphite  Coated  Rocket  Nozzle  Throat  Inserts,"  Final  Report  to 
Atlantic  Research  Corporation,  Cont.  No.  P.  O.  78916,  TRW  Systems  Group, 
Redondo  Beach,  Calif.  , Aug.  28,  1973. 


10 


presents  the  exact  numbers  used  in  generating  the  plots  and  provides  the  material 
needed  to  perform  more  critical  comparisons.  Appendix  A also  contains  tables 
of  % error  for  selected  cases  (Tables  A20-A27).  These  tables  clearly  define 
the  amount  of  absolute  and  relative  error  for  each  of  the  codes. 

In  the  following  discussion  all  displacements  will  be  nodal  point  displace- 
ments (midside  displacements  are  included  in  TEXGAP  results).  On  the  other 
hand,  the  stresses  and  strains  in  all  cases  are  given  for  the  center  of  the  element. 
Since  the  TEXGAP  code  is  the  only  one  which  does  not  output  stresses  and  strains 
at  the  element  center  and  had  to  be  calculated  by  hand  by  averaging  the  output 
values. 

The  locations  of  the  stresses  and  strains  output  from  TEXGAP  depend  upon, 
the  element  type  which  was  used.  For  the  Quad  element  they  are  the  midpoint 
of  the  four  nodal  points.  For  either  element  the  center  element  values  of  stress 
and  strain  were  calculated  by  summing  the  values  at  all  four  points  and  dividing 
by  4.  This  averaged  value  is  the  one  which  appears  in  the  results. 

The  results  for  problems  1 and  2 include  the  exact  solutions.  The  exact 
solutions  were  obtained  from  the  classical  equations  for  pressurized  cylinders 
and  spheres  and  are  listed  at  the  beginning  of  Appendix  A. 

The  data  for  problem  number  1 appear  in  Tables  A!  -A!  and  plots  1-3  in 
Appendix  A.  As  can  be  readily  seen  the  agreement  between  all  codes  and  the 
accuracy  of  each  code  are  very  good  for  this  problem.  It  is  especially  good  for 
radial  displacement  and  strain.  The  TEXGAP  Quad  H element  seems  to  provide 
an  upper  bound  for  all  cases  and,  there  is  a tendency  for  SAAS  III  to  be  lower 
bound,  at  least  in  the  circumferential  direction.  This  is  particularly  apparent 
in  plots  2 and  3. 

The  omission  of  data  for  NONSAP  in  the  results  for  c ircumferential  and 
radial  strains  is  not  an  oversight.  In  its  present  form  the  NONSAP  code  doe6 
not  output  strains. 

Tables  A4-A8  and  plots  4-8  of  Appendix  A contain  the  data  for  problem  2. 

It  was  at  this  point  that  problems  with  the  axisymmetric  solution  in  NONSAP 
were  discovered.  As  a result  NONSAP  was  not  used  for  the  remainder  of  the 
evaluation. 
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Again  the  agreement  between  codes  and  the  accuracy  of  each  code  are  very 
good.  The  tendency  for  TEXGAP  Quad  8 and  SAAS  III  to  form  respectively 
upper  and  lower  data  bounds  is  also  present  in  this  problem  (see  plots  5 and  7). 

An  unexpected  result  observed  in  this  problem  is  the  quite  good  results 
obtained  with  NEPSAP  for  displacements  and  strains  using  a Poisson's  ratio 
of  .4999.  This  element  is  not  reformulated  in  the  sense  that  the  TEXGAP  quad 
element  is  and  yet  the  results  obtained  with  it  rivals  the  TEXGAP  quad  element 
for  accuracy  (see  plots  4 and  6). 

Based  on  the  exact  solutions  for  this  problem,  the  stresses  in  a given 
direction  should  be  the  same  for  any  value  of  Poisson's  ratio.  Only  one  code, 
TEXGAP,  can  be  compared  in  this  manner.  Plot  8 shows  that  the  results 
obtained  from  the  TEXGAP  quad  element  are  very  good  for  this  case. 

The  results  of  problem  number  3 are  presented  in  plots  9-16  and 
Tables  A9-A16.  This  problem  is  common  for  propellant  grains, but  in  this 
case  the  small  L/D  should  be  noted. 

For  problem  3 the  results  for  NEPSAP  do  not  agree  as  well  with  the  other 
codes  for  either  value  of  Poisson's  ratio.  The  difference  becomes  quite 
significant  in  some  cases  (e.  g.  , circumferential  strain,  plots  12  and  13). 

There  are  two  possible  causes  for  the  error  observed,  error  in  the  thermal 
calculations  or  problems  with  the  symmetry  boundry  condition.  Since  the 
results  for  problem  4,  also  a cool  down  problem,  do  not  show  the  same  amount 
of  error  and  considering  that  amount  of  relative  error  increases  with  proximity 
to  the  plane  of  symmetry,  the  boundry  condition  is  most  probably  the  cause. 

There  are  two  schools  of  thought  concerning  strain  output  from  the  codes. 
One  school  of  thought  regards  the  dimensional  change  due  to  free  thermal 
expansion  as  a thermal  strain,  and  the  output  strain  value  includes  this  strain. 
The  other  school  of  thought  considers  the  free  thermal  expansion  to  be  separate 
from  the  strain,  and  therefore  outputs  purely  mechanical  strains.  According  to 
the  first  school  of  thought,  a bar  of  material  subjected  to  a temperature  change 
while  held  between  infinitely  rigid  walls  and  with  no  other  loads  acting  undergoes 
a total  strain  of  zero.  According  to  the  second  school  of  thought,  the  strain  in 
the  bar  is  equal  to  the  free  thermal  expansion  which  would  occur  without  the 
restraint  (but  opposite  in  sign).  Plot  13  shows  that  the  NEPSAP  and  TEXGAP 
Quad  elements  Include  thermal  strain  in  the  output  values. 


16 


For  propellant  grain  problems,  the  stresses  on  a radial  plane  cutting  the 
grain  are  perhaps  more  significant.  A plane  was  passed  through  the  grain  at 
Z = 2.45  inches  and  values  for  all  three  stresses  were  compared  (plots  14,  15, 
and  16).  it  should  be  noted  ag<:in  that  a Poisson's  ratio  of  . 30  was  used  so  that 
all  codes  could  be  romr.»  ••er)  as  can  be  seen,  relative  differences  between  the 
codes  have  increased,  but  the  agreement  for  all  stresses  is  still  very  good. 

Results  for  the  SIC/PG  cool  down  problem  are  shown  on  plots  17,  18,  and 
19,  and  in  Tables  A17-A19.  The  plot  of  circunvferential  stress  (plot  18)  does 
not  contain  data  for  NEPSAP.  This  is  due  to  the  problem  associated  with  the 
addition  of  the  orthotropic  material  model  in  the  axisymmetric  solution.  This 
was  discussed  previously  in  the  section  describing  the  installation  of  the  codes. 

The  results  from  NEPSAP  for  this  problem  would  be  expected  to  show 
similar  relative  errors  to  those  observed  in  the  preceeding  cool  down  problem. 
However,  examination  of  plots  17,  18,  and  19  do  not  show  the  amount  of  difference 
observed  before.  Indeed,  the  results  for  each  of  the  codes  agree  very  well.  Two 
possible  explanations  exist  to  account  for  this  discrepancy.  The  first  is  the 
magnitudes  of  the  numbers.  The  matrix  solution  may  be  having  trouble  with  the 
small  numbers  in  the  previous  grain  cool  down  problem.  Secondly,  this  problem 
does  not  have  a rigid  boundary  condition  restraining  the  outside  diameter. 

Next  to  accuracy,  run  time,  i.  e.  cost,  is  the  most  important  parameter  to 
the  analyst.  Table  3 shows  execution  times  in  seconds  consumed  by  the  central 
processor  (CP  time)  for  each  of  the  codes  in  each  problem.  Also  shown  is  the 
number  of  equations  which  had  to  be  solved  for  each  problem. 

The  fact  which  immediately  stands  out  is  the  amount  of  time  used  by  TEXGAP. 
However,  the  numbers  shown  in  this  table  are  not  completely  fair  to  the  TEXGAP 
code.  The  table  shows  run  times  for  the  same  number  of  elements  that  was  used 
for  the  other  codes.  With  the  elements  in  the  TEXGAP  code,  it  is  possible  to 
use  a much  coarser  mesh  and  still  obtain  results  of  similar  accuracy  obtained 
by  the  remaining  codes  using  the  grids  shown.  These  elements  would  substan- 
tially reduce  the  run  time  for  TEXGAP. 

A second  item  which  can  be  concluded  from  this  table  concerns  NEPSAP.  It 
is  generally  accepted  that  a general  3-D  code  is  more  expensive  to  run  than  the 
2-D  axisymmetric  codes  regardless  of  the  problem  being  solved.  This  is  caused 
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Table  3.  Run  Times  for  Codes 


Plane  Stress  Cylinder 

Number  of  Equations  = 144 
NEPSAP 
NONSAP 
TEXGAP  QUAD 
QUAD  8 

SAAS  III 


Sphere 

Number  of  Equations  = 144 
NEPSAP 
TEXGAP  QUAD 
QUAD  8 

SAAS  III 


Cylinder  Under  Shrinkage 

Number  of  Equations  = 460 
NEPSAP 
TEXGAP  QUAD 
QUAD  8 

SAAS  III 

P.  G.  Cool  Down 

Number  of  Equations  = 525 
NEPSAP 

TEXGAP  QUAD  8 
SAAS  III 


Run  Time  (C.  P.  sec) 


15.  86 
5.  82 
140. 38 
71. 90 
1 1. 05 


14.  51 
139.  58 
71.  75 
1 1. 82 


26.  99 
266. 02 
139. 93 
17.47 


30.  33 
159.  08 
21.43 
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by  the  use  of  sophisticated  solution  techniques  and  core  requirements  needed  for 
higher  order  analyses  (i.  e.  3-D,  nonlinear,  etc.  ) to  solve  problems  which  do 
not  require  such  sophistication.  Although  it  is  true  that  the  run  times  shown  are 
greater  for  NEPSAP  than  for  SAAS  III,  the  difference  is  not  as  great  as  might 
be  expected.  Considering  also  that  NEPSAP  required  only  2000g  locations  more 
to  load  than  SAAS  III,  the  operating  costs  are  very  close. 

USEABILIT Y OF  CODES 

The  remarks  on  the  useability  of  the  codes  is  written  from  the  point  of  view 
that  one  person  is  doing  both  the  job  of  the  analyst  who  writes  up  the  data  input 
in  its  proper  format  and  the  keypuncher  who  must  punch  the  cards. 

The  simplest  of  the  four  codes  from  a user  standpoint  is,  without  question, 
TEXGAP.  Its  free-field  input  makes  data  input  much  easier  than  any  of  the  other 
codes.  The  only  restrictions  which  must  be  followed  are  the  order  of  the  data 
on  each  record  (card)  and  to  separate  each  item  on  the  record  with  a comma. 

The  mesh  generation,  element  and  boundary  condition  assignments  are  also 
much  simpler  than  the  other  three  codes  and  possibly  simpler  than  any  other 
code  available  to  the  analyst. 

The  biggest  problem  in  using  TEXGAP  occurs  when  inputting  material 
properties  for  orthotropic  materials.  TEXGAP  requires  manual  input  of  the 
actual  stiffness  constants  for  this  case  instead  of  Young's  Moduli  and  Poisson's 
ratios  for  the  material  directions. 

SAAS  III  ranks  second  in  useability.  Its  method  for  data  input,  which  is 
typical  for  most  codes,  requires  the  user  to  use  specific  columns  for  the  data. 

The  mesh  (if  generated  by  the  code)  is  specified  by  supplying  line  segments, 
straight  or  curved,  which  bound  the  region  being  modeled.  The  elements  can 
only  be  generated  internally  in  lines  parallel  to  lines  of  constant  I or  J values. 

This  requires  several  cards  of  input  data  to  generate  even  the  simplest  geometries. 
Computation  of  the  mesh  in  desired  areas  is  possible,  but  the  procedure  is  not 
simple  to  understand.  It  would  require  substantial  work  to  become  familiar 
with  its  use  and  probably  some  trial  and  error  even  then  to  get  the  desired  mesh. 

Several  questions  concerning  some  of  the  job  control  parameters  came  up 
when  learning  to  use  SAAS  III.  For  example,  it  is  not  clear  to  the  user  that  a 
0 (zero)  for  the  plane  stress/plane  strain  option  gives  an  axisymmetric  solution. 
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Similarly  for  the  stop  parameter,  it  isn't  obvious  that  a 0 (zero)  makes  the 
program  do  a complete  solution.  These  things  are  minor  but  frustrating  to  the 
user.  Much  time  digging  through  the  program  listing  could  be  avoided  by 
providing  a complete  set  of  options,  including  defaults,  for  all  input  parameter 
in  the  input  guide.  One  final  item  which  caused  some  question  was  the  magnitudes 
of  strain  which  the  program  outputs.  After  searching  through  the  program 
listing,  it  was  determined  that  the  strains  are  output  in  percent;  but  neither  the 
user  manual  nor  the  output  states  that  fact.  The  user's  guide  is  well  written 
and  is  reasonably  easy  to  understand.  The  only  difficulties  encountered  in  its 
use  were  associated  with  the  loads  data.  There  are  several  parameters  on 
different  data  cards  which  control  the  loads  data  and  the  instructions  lack  a little 
detail  concerning  their  use.  Also  it  was  discovered  that  two  blank  cards  are 
needed  at  the  end  of  the  data  to  get  a normal  termination  of  the  run. 

NONSAP  was  used  very  little  in  this  evaluation  due  to  the  problems  encoun- 
tered with  the  axisymmetric  solution.  Two  items  were  discovered  concerning 
the  use  of  NONSAP.  The  first  item  is  that  any  load,  thermal  or  mechanical, 
must  be  converted  to  concentrated  nodal  loads  externally.  This  considerably 
complicates  the  input  of  loads.  The  second  item  is  the  lack  of  strain  output. 

DEFICIENCIES  OF  CODES 

In  order  to  be  able  to  describe  the  deficiencies  of  the  current  codes,  the 
needed  capabilities  should  first  be  defined.  The  features  desired  in  a general 
purpose  nozzle  code  were  discussed  in  detail  in  Volume  I of  this  report.  To 
summarize  the  desired  code,  it  would  be  a finite  element  code  with  2-D  axisym- 
metric and  3-D  elements.  It  must  be  capable  of  analyzing  a variety  of  materials 
that  range  from  isotropic  to  fully  anisotropic  in  character  .id  allow  temperature 
dependent  properties.  It  must  be  capable  of  handling  axisymmetric  as  well  as 
non-axisymmetric  loads.  These  loads  would  include  concentrated  pressure  and 
thermal  loads  plus  initial  displacements  and  strains.  The  code  must  also  allow 
material  and  kinematic  (geometric)  nonlinearities.  And  finally,  it  must  be 
nonproprietary.  These  capabilities  would  form  an  adequate  nozzle  code. 

There  are  several  other  capabilities  deemed  very  desirable  for  this  code, 
although  they  are  not  absolutely  necessary.  High  on  this  list  would  be  generation 
schemes  for  nodes  and  elements  plus  graphics,  either  on  line  or  paper  plots, 
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for  developing  the  grids  and  displaying  output.  A system  that  allows  multiple 
I - J grids  and  makes  it  possible  to  alter  the  grid  of  separate  components  without 
disturbing  the  entire  model  would  be  useful.  Boundary  conditions  which  allow 
adjacent  components  to  slip  or  model  gap  closure  would  be  very  desirable.  Also 
such  things  as  compatible  high  and  low  order  elements  to  allow  maximum 
efficiency  in  gridding,  in  dynamic  analysis  capability,  and  in  companion  thermal 
analyses  would  be  useful  items. 

Speaking  in  general  terms,  the  chief  deficiency  for  all  codes  appears  to  be 
the  lack  of  a nonlinear  anisotropic  material  model.  The  DOASIS  code  is  the  only 
one  of  the  codes  studied  which  has  a theoretically  sound  nonlinear,  small  strain, 
anisotropic  model.  This  model  assumes  the  material  has  three  mutually  orthogonal 
principal  axes,  i.  e.  the  material  must  be  orthotropic.  The  small  strain  assump- 
tion is  probably  sufficie  t for  most  of  the  state  of  the  art  materials  since  most 
have  a low  strain  to  failure. 

None  of  the  codes  currently  used  allow  sliding  contact  or  gap  closure  bound ry 
conditions.  These  would  be  particularly  useful  for  modeling  throats  made  of  PG 
washers,  but  there  are  numerous  places  in  present  nozzle  designs  which  display 
these  types  of  behavior. 

Of  the  four  codes  which  were  run,  NEPSAP-  has  the  fewest  deficiencies,  but 
it  still  lacks  some  important  capabilities.  NEPSAP  does  not  have  an  anisotropic 
axisymmetric  element  and  also  lacks  a nonlinear  anisotropic  model. 

As  for  NONSAP,  the  deficiencies  are  also  in  the  lack  of  anisotropic  material 
models.  Most  of  the  nonlinear  capabilities  in  NONSAP  are  for  isotropic  materials 
and  are  useable  for  2-D  axisymmetric  elements  only.  In  addition,  none  of  the 
material  models  allow  temperature  dependent  properties.  Finally,  the  code 
permits  only  concentrated  model  loads  which,  in  the  case  of  2-D  axisymmetric 
analysis,  must  be  axisymmetric. 

SAAS  III  is  a 2-D  axisymmetric  code  and  therefore  has  no  3-D  capabilities. 

The  code's  nonlinear  capability  consists  solely  of  a bilinear  curve  fit  on  Young's 
modulus.  This  model  must  be  used  with  much  care  given  to  determining  its 


*This  refers  to  the  version  of  NEPSAP  obtained  from  ASIAC. 
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applicability.  The  level  of  anisotropy  extends  only  to  orthotropic  materials  and 
they  must  be  axisymmetric.  Finally,  SAAS  III  does  not  allow  non-axisymmetric 
loadings. 

Like  SAAS  III,  TEXGAP  is  a 2-D  axisymmetric  code  and  does  not  handle 
3-D  analysis,  but  it  has  no  nonlinear  capability  at  all.  Again,  like  SAAS  III, 
its  level  of  anisotropy  extends  only  to  orthotropic  materials  and  they  must  be 
axisymmetric.  It  allows  concentrated  loads  only  by  approximation  through  a 
Fourier  series  expansion  technique.  Of  most  significance  TEXGAP  does  not 
allow  temperature  dependent  material  properties  nor  analysis  of  a problem 
containing  thermal  gradients. 


22 


SECTION  IV 


CONCLUSIONS 

The  only  surprising  result  uncovered  in  the  program  is  the  reasonably  good 
values  for  displacements  and  strains  obtained  by  NEPSAP  with  a Poisson's  ratio 
near  one-half.  The  accuracy  of  the  results  obtained  by  each  of  the  four  codes  on 
all  of  the  problems  considered  in  this  project  is  good.  This  is  true  for  the 
relative  agreement  between  the  codes  and  for  the  absolute  accuracy  of  each  code 
in  the  problems  with  exact  solutions. 

None  of  the  codes  contained  all  nor  most  of  the  features  desirable  for  a 
general  purpose  nozzle  code,  but  many  of  these  features  were  addressed  by  at 
least  one  of  the  codes  if  only  for  the  isotropic,  elastic  case.  This  represents 
a basic  understanding  by  the  code  developers  for  the  needs  of  the  analysts  in  the 
nozzle  industry  and  provides  a framework  on  which  to  build  better  codes. 

The  calculated  strain  values  from  any  computer  code  must  be  viewed  cau- 
tiously when  solving  problems  involvong  thermal  expansion.  The  output  values 
of  strain  may  or  may  not  contain  the  thermal  strain  depending  on  which  code  is 
being  used.  In  one  case  using  the  TEXGAP  code,  these  values  even  depended 
on  which  element  was  being  used. 

Based  on  the  results  of  this  program  it  is  not  necessarily  true  that  the  higher 
order  elements  give  "better"  results.  While  it  is  true  that  the  higher  order 
elements  used  in  a coarser  grid  will  provide  comparable  results  with  the  standard 
elements,  a blanket  statement  about  better  results  can't  be  made.  The  deter- 
mination of  the  point  at  which  these  elements  begin  to  make  visable  payoffs  when 
considering  accuracy  versus  grid  coarseness  and  run  times  deserves  more 
investigation. 

One  final  item,  and  one  that  is  generally  accepted  as  true  but  did  not  prove 
to  be  very  significant,  is  the  increased  cost  which  must  be  paid  to  run  a general 
3-D  code,  over  that  of  the  standard  2-D  axisymmetric  codes,  to  analyze  2-D 
problems.  If  NEPSAP  is  accepted  as  a general  3-D  code  and  the  problems 
analyzed  in  this  project  considered  as  reasonably  similar  to  typical  problems, 
then  the  concept  of  user  cost  needs  to  be  re-evaluated. 
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SECTION  V 


APPENDIX  A 


The  results  of  the  first  two  problems,  internally  pressurized  cylinder  and 
sphere  respectively,  were  compared  to  results  obtained  from  the  classical 
solutions  for  these  two  problems.  For  convenience  the  equations  used  to  deter- 
mine the  exact  solutions  are  given  below.* 

Internally  Pressurized  Cylinder 
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Theory  of  Elasticity,  Third  Edition,  Timoshenko  and  Goodier,  pp.  70-71  and 
pp.  392-395. 
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Internally  Pressurized  Sphere 
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Table  A-4 

Sphere  Under  Internal  Pressure  (Nu  = .4999) 
Radial  Stress  (P.  S.  I.  ) Vs  Radial  Position  (inches) 


Radial 

Exact 

TEXGAP 

Position 

Solution 

Quad  Flem. 

1.  1 

-74. 86 

-74. 93 

1. 35 

-39. 99 

-40. 49 

1.7 

-19.47 

-19.  79 

2.  15 

-9.  06 

-9.  27 

2.  7 

-4.  03 

-4.  09 

3.  35 

-1.  58 

-1. 60 

4.  1 

-. 3575 

-. 3795 

-6 


Table  A-fe 

Sphere  Under  Internal  Pressure  (Nu  = .4999) 
Circumferential  Strain  (in.  /in.  ) Vs  Radial  Position  (inches) 


Radial 

Position 

Exact 

Solution 

NEPSAP 

T EXGAP 
Quad.  Elem. 

1.  1 

. 05697 

. 05607 

. 05626 

1.  35 

. 03082 

. 03054 

. 03088 

1.7 

. 01543 

. 01528 

.01554 

2.  15 

. 007630 

. 007527 

. 007672 

2.  7 

. 003853 

. 003774 

. 003859 

3.  35 

. 002017 

. 001961 

. 002011 

4.  1 

. 001100 

. 001067 

. 001092 

-8 
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Table  A-8 

Sphere  Under  Internal  Pressure 
Radial  Stress  (P.  S.  I.  ) Vs  Radial  Position  (inches) 


Radial 

Position 


1.  1 


1.  55 


1.7 
2.  1 5 


2.  7 


5.  55 


4.  1 


Exact 

Solution 


-74.  86 


-39.  99 


-19.47 


-9.  06 
-4.  03 
-1.  58 


Nu  .4999 
TEX GAP 
Quad.  Elem. 


-74. 93 


-40. 49 
-19. 79 


Nu  = .30 
TEX GAP 
Quad.  Elem. 


-75.  1 1 


-40.  63 


19.  88 


3575 


Table  A-9 


Cylinder  Under  Shrinkage  (Nu  = .4999) 

Radial  Displacement  (inches)  Vs  Axial  Position  (inches) 


Table  A-  10 


Cylinder  Under  Shrinkage  (Nu  = . 30) 

Radial  Displacement  (inches)  Vs  Axial  Position  (inches) 


Axial 
Pos  ition 

NFPSAP 

0. 

. 008267 

. OS 

TEXGAP 

TEXGAP 

Quad  8 Elem. 

Quad  Elem. 

. 007833 

. 007826 

. 008318 

. 008310 

. 008804 

. 008795 

. 009535 

. 009527 

. 01027 

.01026 

.01125 

.01124 

SAAS  III 
. 007738 

. 008709 

. 01018 


.01456 

.01455 

.0158b 

. 01585 

.01707 

.01707 

.01836 

. 01835 

. 01950 

.01949 

. 02063 

. 02062 

. 02158 

. 02158 

, 

. 02246 

. 02245 

.02315 

. 02315 

. 02375 

. 02375 

. 024  1 9 

. 02419 

2 
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Table  A-12 

Cylinder  Under  Shrinkage  (Nu  = .4999) 
Circumferential  Strain  (in.  /in.  ) Vs  Axial  Position  (inches) 


Axial 

Position 

NEPSAP 

T EXGAP 
Quad.  Elem. 

Thiokol 

Isopara.  Elem. 

. 05 

. 01432 

.01112 

.01144 

. 175 

. 01731 

.01443 

.01470 

. 35 

. 02143 

.01910 

. 01935 

. 575 

. 02694 

. 02508 

. 02532 

. 85 

. 03284 

. 03205 

. 03230 

1.  175 

. 03923 

. 03948 

. 0 3974 

1.  55 

. 04460 

. 04652 

. 04694 

1. 975 

. 04923 

.05299 

. 05327 

2.45 

. 05207 

.05799 

. 05826 

2.  975 

~ .05374 

.06145 

. 06169 

3.  55 

.05338 

. 06324 

. 06346 
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03270  . 02957  .03274 
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19.  51 


Table 


Table  A- 17 

SiC/PG  Cooldown  (Throat  plane) 
Displacement  (inches)  Vs  Radial  Position 


NEPSAP 


-.005401 


TEXGAP 
Quad  8 Flem. 


005395 


005447 


005499 


005600 


005700 


005800 


005900 


006000 


006100 


-.006199 


-. 006299 


-.006399 


-. 006498 


-. 006548 


-. 006598 


-.006635 


-. 006672 


-. 006764 


-. 006856 


-. 007218 


-. 007581 


-. 007947 


-. 008314 


-. 008683 


-.009053 


-. 009424 


-. 009796 


-.009833 


-. 009870 


(inches) 


SAAS  III 


-. 005400 


-.005501 


-. 005701 


-. 005901 


006101 


-.  006301 


-.006500 


-. 006600 


-. 006673 


-. 006857 


-.007629 


-. 008367 


-. 009086 


-. 009798 


-. 009875 


Table  A-  18 

SiC/PG  Cooldown  (Throat  plane) 
Circumferential  Stress  (P.  S.I.  ) Vs  Radial  Position  (inches) 


A-20 


a 


Table  A-19 

SiC/PG  Cooldown  (Throat  Plane) 
Circumferential  Strain  (in. /in.  ) Vs  Radial  Position  (inches) 


Radial 

Position 


. 8165 


. 8355 


. 8605 


. 8855 


5 


. 9355 


1 . 4 30 


1.494 


NEPSAP 
Thermal  Strain 
Removed 


. 000420 


T EXGAP 
Quad.  8 Elem. 


SAAS  III 


. 000227 


-.  000669 


-.  00040  3 


-.000385 


. 0004308 

. 00042 

. 0004001 

. 00039 

. 0003614 

. 00035 

.0003256 

. 00032 

. 0002923 

. 00029 

. 0002609 

. 00025 

. 0002384 

. 00023 

-. 0006628 

-. 00067 

-. 0006432 

-.  00065 

-. 0005844 

-. 00058 

0005063 

-. 00050 

-. 0004480 

00045 

0004034 

-. 00040 

-. 0003823 

i 

-. 00039 

-2 
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Table  A-24 

Relative  Error  (%)  Based  on  TEXGAP  Quad  Element 
Cylinder  Under  Shrinkage  (Nu  = . 30)  - Radial  Displacement 


Axial 

Location 

NEPSAP 

TEXGAP 
Quad  8 Elem. 

SAAS  III 

0. 

5.  6 

. 1 

8 

. 7 

1.  3 

. 1 

3 

2.  2 

1.8 

0. 

. 2 

3.  85 

1.  7 

0. 

0. 

The  TEXGAP  Quad  Element  was  chosen  as  the  norm  since  it 
consistently  had  the  smallest  error  in  the  first  two  problems. 
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Table  A-  25 

Relative  Error  (%)  Based  on  TEXGAP  Quad  Element 
Cylinder  Under  Shrinkage  (Nu  = .30)  - Circumferential  Stress 


Rad  ial 
Location 

NEPSAP 

TEXGAP 
Quad  8 Elem. 

SAAS  III 

1. 05 

-7.  0 

. 2 

-.  7 

2.  25 

-1.6 

. 7 

1.4 

3.  95 

-.  5 

0. 

. 3 
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Table  A-26 


Relative 
SiC/PG  Cooldown 


Error  (%)  Based  on  SAAS  III* 

(Throat  plane)  - Radial  Displacement 


Rad  ial 
Location 

NEPSAP 

TEXGAP 
Quad  8 Elem. 

.810 

0. 

-.  1 

. 848 

. 1 

0. 

.923 

. 1 

0. 

. 960 

0. 

0. 

1.004 

0. 

0. 

1.  372 

0. 

-.  4 

1. 500 

0. 

0. 

The  TEXGAP  Quad  element  could  not  be  used  for  this  problem  due 
to  the  lack  of  an  Orthotropic  material  model.  Therefore,  SAAS  III 
was  chosen  as  the  basis. 
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Figure  A-5.  Ci  reunite  rent  iai  Stress  Vs  Radial  Location 
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Figure  A-12.  Circumferential  Strain  Vs  Axial  Locat 


Figure  A-13.  Circumferential  Strain  Vs  Axial  Location 
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